Threshold for everlasting initial memory in equilibration processes 
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Conventional wisdom indicates that initial memory should decay away exponentially in time for 
general (noncritial) equilibration processes. In particular, time-integrated quantities such as heat 
are presumed to lose initial memory in a sufficiently long-time limit. However, we show that the 
large deviation function of time-integrated quantities may exhibit initial memory effect even in the 
infinite-time limit, if the system is initially prepared sufficiently far away from equilibrium. For a 
Brownian particle dynamics, as an example, we found a sharp finite threshold rigorously, beyond 
which the corresponding large deviation function contains everlasting initial memory. The physical 
origin for this phenomenon is explored with an intuitive argument and also from a toy model analysis. 
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Hot coffee gets colder and iced coffee gets warmer at 
room temperatures. These phenomena are the examples 
of equilibration processes and can be generalized as the 
following situation; a system with initial temperature T s 
is in thermal contact with a heat bath with temperature 
T;,. Then, the system gradually deviates from its ini- 
tial state and approaches to the final equilibrium (EQ) 
state which is determined by the heat bath. Here, the 
initial distance from final equilibrium is parameterized 
by the temperature ratio (3 = T^/T s . The relaxation 
process is usually exponentially fast, so the memory of 
the initial temperature will be lost for average values of 
most physical observablcs after a characteristic relaxation 
time. However, the initial memory can often survive in 
the tail part (rare-event region) of a probability distribu- 
tion function (PDF) even in the long-time limit. 

What about time-integrated quantities such as heat, 
work, or entropy production, which are the key quanti- 
ties for nonequilibrium (NEQ) fluctuation theorems [H- 
5,]? These accumulated quantities are also affected by a fi- 
nite transient period, but their average values increase (or 
decrease) linearly in time asymptotically in NEQ steady 
state. Therefore, in a sufficiently long-time limit, our 
conventional wisdom may lead us to expect that they 
will lose all initial memory (independent of (3). Never- 
theless, in this Letter, we show rigorously that this is false 
wisdom for time-integrated quantities and, in particular, 
corresponding large deviation functions depend strongly 
on the initial condition (/3) even in the infinite-time limit. 
More surprisingly, there exists a sharp threshold for /3 -1 
in general, only beyond which the initial memory lasts 
forever. 

In literatures, there have been some reports that ini- 
tial conditions can affect the large deviation function in 
the long time limit [B^l • For example, van Zon and 
Cohen [6] showed that heat transfer Q in a driven har- 
monic oscillator in contact with a heat bath violates the 
fluctuation theorem even in the long-time limit, starting 
initially from EQ. In contrast to work, heat is known to 



satisfy the fluctuation theorem, only starting with a uni- 
form distribution (infinite-temperature initial state) [Io| . 
Thus, their finding can be interpreted as an everlasting 
initial memory effect in the large deviation function for 
heat. 

In this Letter, we consider heat transfer during the 
equilibration process of a simple Brownian particle and 
investigate initial memory effects systematically in the 
long-time limit. The Brownian particle dynamics is de- 
scribed by the Langevin equation 



(1) 



where v is the velocity of the particle, 7 is the dissi- 
pative coefficient, and £ denotes a random white noise 
satisfying (£(t)£(t')) = 2DS(t — r'). Here, we set the 
particle mass m = 1 for convenience and the heat bath 
temperature T\, = D/j. Initially, the system is prepared 
in EQ state with the Boltzmann distribution at tempera- 
ture T s = Tb/ [3. And then, the thermal contact is formed 
at time r = between the system and the heat bath, and 
maintained until final time r = t. 

Time-integrated heat flow between the system and the 
heat bath can be decomposed into the dissipated energy 
flow Qd from the system into the heat bath and the in- 
jected energy flow Qi in the other way around Q: 



dr jv and 



dr £,v. 



(2) 



Even if the system reaches EQ in the long-time limit, 
each of (Qd) and (Qi) increases linearly in time t indef- 
initely with their difference representing the system en- 
ergy change (AE) = ^[(v 2 (t)) - (v 2 (0))}, which is finite 
for nonzero (3. As expected, there will be no net heat 
flow at EQ. 

We first study the PDF of the (average) dissipated 
power, Ed = Qd/t, and later the injected power, Si = 
Qi/t. To calculate the PDF, P(ed), it is convenient to 
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consider its generating function 

/oo 
de d P(e d )e- xt ^, (3) 
-oo 

which is the Fourier transform of P(s d ). The generating 
function can be calculated exactly by the standard path 
integral method 0, EH • With the initial Boltzmann dis- 
tribution Pinit(f(0)) ~ exp[— u 2 (0)/(2T s )] at temperature 
T s = D/(j/3), we find 



7Td 



(A) = e*/ 2 f coshmt + LLhll si 



-1/2 



sinh 7777; 



(4) 



with dimcnsionless parameters A = 2D A/7 and 7/ = 

\/r+2A. 

The inverse Fourier transform of Eq. Q yields the 
PDF in terms oi £ d = £ d / D as 
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exp 



f (EdA + 1) 



cosh rj^jt + sinh 7777; 



(5) 



For large t, the above integration can be carried out by 
the saddle point approximation. However, care should 
be taken due to the presence of the branch cut. Here, we 
take the branch cut on the rcal-A axis where 



/ = cosh ijjt 



1 + A//3 



sinh 777/; 



(0) 



becomes negative, see Fig. [TJ 

We locate the branch points for large t, which depend 
on /3. Note that 77 is real and positive for A > — i, while 
r] becomes pure imaginary for A < — |. For /3 > |, / > 
and we have no branch points for A > — i. Instead, we 
find them in the region of A < — | and the largest one is 
denoted by A^ ~ — i + 0(t~ 2 ). In the t — > 00 limit, the 
branch point A^ = — i has no /3 dependence. For /3 < ^ , 
in contrast, / can become negative for A > — i and we 
find the branch point approaching A^ = —2/3(1 — /3) as 
t — > 00. Locations of A^'s and branch cuts are shown 
in Figs. [l{a) and (b). It turns out that the branch-cut 
structure plays a crucial role in determining the everlast- 
ing initial memory effect. 

From Eq. ([5]), one may easily expect for large t 



P(e d ) ^ exp [th t (e d )] 



(7) 



with the large deviation function (LDF) h(ed) = 
lim t _ i . 00 h t (sd)- We first calculate the LDF using the 
saddle point method in the presence of the branch-cut 
structure found as above. For large t, Eq. becomes 
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FIG. 1. (Color online) (a) and (b) show the branch-cut struc- 
ture of TTd{\) on the complex A plane for j3 > 1/2 and /3 < 1/2, 
respectively, (c) and (d) show the branch-cut structure of 
7Ti(A) for P > 1/4 and /3 < 1/4, respectively. Wiggled lines 
denote branch cuts. A head of an each arrow locates at its 
asymptotic value of the respective branch point as t — > 00. 



/ ~ \e^ 1 {\ + ■ The saddle point X* is given by 

the solution of the following equation: 



d_ 

<fx 



= 0, (8) 



where the logarithmic term is included because it may 
become very large in the vicinity of A = —2/3(1 — /3) for 
/3<i. 

For /3 > i, we find a solution (saddle point) on the 
real-A axis which is outside of the branch cut as 



1 



A* 



as t — > 00. Then, the LDF for (3 > 1/2 becomes 



M^) = + 1 -v*) 



-is, 



(9) 



(10) 



where if = y 1 + 2A^ and the logarithmic term is negli- 
gible. As P(sd) = for Ed < 0, the LDF is defined only 
for e~d > 0. This LDF has no /3 dependence but is deter- 
mined only by the heat bath properties (7, D). Thus, we 
call Eq. (fTOj) the heat-bath characteristic curve (HBCC). 

For /3 < i, the saddle point location exhibits a non- 
analytic behavior as function of Ed, due to the interplay 
of the saddle point and the branch point. When e d < 
(1 — 2/3) -1 , the saddle point AJ given by Eq. (0) is located 
to the right of the branch point, X* d > X d = -2/3(1 — /3), 
in the t — s- 00 limit. Thus, the LDF h(£ d ) is identical 
to the HBCC in Eq. When e d > (1 - 2/3) -1 , the 

saddle point approaches the branch point asymptotically 
from the right side, due to the divergence of the loga- 
rithmic term in Eq. at the branch point. However, 
as this approach is not exponentially fast in time, the 
dominant contribution to the LDF comes from the con- 
ventional first term in Eq. ([8|) at the asymptotic saddle 



point A* = A d 



-2/3(1 — /3). Summarizing for /3 < i, 
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the LDF is 

h(e d ) 



7/3 [(1 - P)e d - 1] 



£d < 1-2/3 
£d > T=2fi 



(11) 



Note that the LDF for large Ed is deformed from the 
HBCC and has the initial condition (/3) dependence, see 
Figs, da) and (b). 

Our results show that, for sufficiently high initial tem- 
peratures (/3 < the initial memory survives for- 
ever in the large Ed region of the LDF and completely 
vanishes below the threshold of the initial temperature 
(P > |). Large dissipated energy Qd is generated by 
the decay of highly energetic particles with energy ~ Qd- 
There are two distinct sources for highly energetic par- 
ticles; (a) heat bath and (b) initial Boltzmann distribu- 
tion, which compete each other. We estimate the prob- 
ability P a and Pb to find a particle to dissipate energy 
Qd from each source, respectively. From the HBCC in 
Eqs. © and (HOD, we find P a ~ cxp[-Q d /(4T 6 )] for 
large E d - O n the other hand, we assume that a parti- 
cle with high initial energy decays by the deterministic 
dynamics of v = —jv. In this case, the dissipated en- 
ergy is Qd = 5^ 2 (0) in the long-time limit. Thus, we 
estimate P& ~ exp[— Qd/T s ] = exp[— PQd/T ] from the 
initial Boltzmann distribution at temperature T s . As a 
result, the HBCC P a dominates over for p > ^ or 
the initial memory dominates, otherwise. As Qd is over- 
estimated in the latter case, the threshold value \ only 
sets its lower bound, which is consistent with the correct 
value P c d = \- 

We also calculate the leading finite-time correction of 
ht(sd) in Eq. (JT]). As the leading correction is 0(lnt/t), it 
yields a power-law type prefactor to the exponential form 
of the PDF, P(sd)- For /3 < g, it is tricky to calculate 
this correction because the saddle point is very close to 
the branch point. In fact, it cannot be obtained through 
a conventional Gaussian integral. Here, we just report 
our result without presenting details [l2| for (3 < i 
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(B) 
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J12) 

where there are three regions: (A) (1 — 2/3) 1 — Sd 
(TtJ-Va; ( B ) ^-(l-2/3)- 1 |«( 7 t)- 1 /2 ; (C)e d -(1- 
2/3) 1 ^s> ("ft) 1//2 . Three constants are given as Cd(P) = 
2^(l-2/3) 7 /4r(I)/(47ryr^), and 
For /3 > |, the PDF 
is given by the same one in (A) of Eq. (|T2|) for all Id > 0. 
The prefactors depend on the initial condition (/3) for 
all cases, as expected, but their power-law exponent in 
terms of £ d changes abruptly from —2 to — | as e d in- 
creases. It is interesting to note that this exponent 
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FIG. 2. (Color online) (a) and (b) are the LDF's of the dis- 
sipated power for /3 > 1/2 and /3 < 1/2, respectively. The 
solid line is the HBCC (see Eq. (fTCT)) '). In (a) x, *, o, and ■ 
are numerical data for ft — 3/4, 1, 2, and 4, respectively, (c) 
and (d) are the LDF's of the injected power for /3 > 1/4 and 
P < 1/4, respectively. In (c) x, *, o, and ■ are numerical 
data for /3 = 1/2, 1, 2, and 4, respectively. Each dashed line 
denotes the analytic line for each /3. All numerical results are 
obtained at t = 100. 



change is very similar to what was found dynamically 
for the PDF of nonequilibrium work in simple linear dif- 
fusion systems [l3j . 

Now, we turn to the injected power, £j = Qi/t. The 
calculation method is almost the same as before. We 
obtain the generating function of the injected power as 



7T<(A) 



o7«/2 



cosh rpft 



1 + A- A 2 /2/3 



'/ 



\ -1/2 

sinh?/7t (13) 



Compared to Eq. ((4|), there is only a parametric differ- 
ence in the coefficient of the hyperbolic sine term. We 
can obtain the PDF of the dimensionless injected power, 
Pili) with £i = Ei/D, by performing the inverse Fourier 
transform of 7Tj(A). 

Similar to the case of the dissipated power, the branch 
points are determined by the equation 



= cosh 777/ + 



1 + A- A 2 /2/3 



sinh?77<. (14) 



We find two relevant solutions of Eq. (|14p ; one is on the 
positive real axis, A^~, and the other is on the negative 
real axis, A~, as shown in Figs.QJc) and (d), respectively. 
In the t — > 00 limit, one can show that A^~ = 2(/3 + y/]3) 
for all /3, while \~ = -\ for p > \ and A t ~ = -2( V/3-/3) 
for p<\. 

By defining the LDF, h(e~i), for ei in the t — > 00 limit 



4 



and through a similar algebra, we find for /3 > 



h{ei 



[1 - (1 + VP)Ei] , Ei < 

I) 2 , > 



4ei V fc 



(15) 



Note that /i(e~i) is defined for all £j. The negative tail 
of P(?i) is affected by the initial condition and the non- 
analyticity of h(si) is present even in the EQ process at 
(3=1. For /3 < ^ i the LDF becomes 



/l(£i) 



- 7 V£[1-(1 
-&(^-l) 2 > 



i 



1+2^ 

[(1 - V?)ei - 1] , 



1+2V3 
< < t^-ts(16) 



e 4 > 



1-2^ 
l-2y^ 



Our results read that the negative tail always depends 
on the initial condition, but the positive tail shows the 
threshold at ftf = j where the initial condition depen- 
dence starts to appear. Note that the threshold value 
varies with the quantity interested. The leading finite- 
time correction is rather complicated, which will appear 
elsewhere 
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To confirm our analytic calculations in Eqs. (|10[l . ([TT]) . 
(|15p . and (|16j) . we performed numerical integrations of 
the Langevin equation, Eq. (TTJ. Here, we set 7 = D = 1, 
and integration time interval At = 10~ 3 . Figure Uta) 
displays numerical data for h(ed) at t = 100 for various 
values of j3 > 1/2. Regardless of /3, all numerical results 
collapse well onto the HBCC as expected from Eq. (fTUl) . 
Slight deviation from the analytic HBCC comes from the 
finite-time effect. We confirmed that the LDF with lead- 
ing finite-time correction (see Eq. (JUJ)) perfectly agrees 
with the numerical data at t = 100 (not shown here). 
Figure [UJb) shows h(ed) for /3 < 1/2. Numerical results 
also agree well with our analytic results in Eq. (flT) . Fig- 
ure l^c) and (d) show the LDF of the injected power for 
P > 1/4 and /? < 1/4, respectively In Fig.^c) the LDF 
for Ei < (1 + 2y/]3)~ 1 does not appear simply because the 
region is outside of the plot range. Meanwhile, the three 
regions are clearly seen in Fig. [U[d), as expected from 
Eq. (HU). 

To understand better the origin of the threshold (3° di- 
viding different phases, we introduce a simple toy model. 
In order to examine the correlation between the initial en- 
ergy and the average power, we define a function E t (i) 
which is the average initial energy of a particle whose 
average (dissipated or injected) power until time t is 
given by e. It is convenient to use the normalized func- 
tion e(e) = E(e)/(E)q, with (E)o the average initial en- 
ergy without any constraint on its power. In Fig. EHa), 
e(ffd) approaches the constant 1 for large t, which im- 
plies no correlation between the initial energy and the 
corresponding dissipation power. Thus, there will be no 
initial condition dependence on the PDF for large Ed- 
In contrast, Fig. [21b) shows divergence of e(£d) m time, 
which indicates that high initial energy is responsible for 




FIG. 3. (Color online) (a) and (b) show the correlation in 
time between the initial energy and the dissipated power for 
j3 > 1/2 and j3 < 1/2, respectively, (c) and (d) show the 
correlation between the initial energy and the injected power 
at time r for /3 > 1/4 and p < 1/4, respectively. A, 0, and o 
are data for r = 1, 10, and 100, respectively. 



large Ed- Thus, the tail of the PDF, P(£d), should be 
dominated by a particle with high initial energy which 
is generated by the initial high-temperature distribution. 
This causes the deviation of h{£d) from the HBCC for 
Ed > (1 — 2/?) _1 when /3 < | (higher initial tempera- 
tures). FigureEJc) shows e(£i) approaches 1 for j3 > 1/4, 
while Figure [3jd) shows its divergence for /? < 1/4. 

In summary, we consider the equilibration process of 
a Brownian particle system, staring from various initial 
temperatures different from the heat bath temperature. 
We calculate the LDF of time-integrated quantities like 
the dissipated energy and the injected energy due to the 
heat bath. Remarkably, we find a finite threshold for the 
initial temperature, only beyond which the LDF contains 
everlasting initial memory. We argue that this is due to 
the competition of highly energetic particles originated 
from the heat bath and from the initial distribution. 
Our simple toy model analysis supports this argument by 
showing that large dissipated energy is generated domi- 
nantly by particles with high initial energy, rather than 
by highly energetic particles randomly generated by the 
heat bath, when the initial temperature is sufficiently 
high enough with respect to the heat bath temperature. 
We expect that our results are applicable to general equi- 
libration or noncquilibrium processes, which implies that 
the rare-event measurements for time-integrated quan- 
tities should be carefully carried out because the initial 
memory may survive forever. 
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